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Abstract. Let {Wi}i£i be a (redundant) sequence of subspaces each being endowed with a 
weight Vi, and let Ti be the closed linear span of the Wi's, a composite Hilbert space. Provided that 
{{Wi,Vi)}i£i satisfies a certain property which controls the weighted overlaps of the subspaces, it is 
called a fusion frame. These systems contain conventional frames as a special case, however they go 
far "beyond frame theory". In case each subspace Wi is equipped with a frame system {fij}j£j by 
which it is spanned, we refer to {{Wi,Vi, {/ij }jg )}ig/ as a fusion frame system. 

In this paper, we describe a weighted and distributed processing procedure that fuse together 
information in all subspaces Wi of a fusion frame system to obtain the global information in Ti. The 
weighted and distributed processing technique described in fusion frames is not only a natural fit 
in distributed processing systems such as sensor networks, but also an efficient scheme for parallel 
processing of very large frame systems. We further provide an extensive study of the robustness of 
fusion frame systems. 
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1. Introduction. Frames, which are systems that provide robust, stable and 
usually non-unique representations of vectors, have been a focus of study in the last 
two decades in applications where redundancy plays a vital and useful role, e.g., filter 
bank theory [HI, sigma-delta quantization signal and image processing 0, and 
wireless communications |19| . 

However, a number of new applications have emerged where the set-up can hardly 
be modeled naturally by one single frame system. They generally share a common 
property that requires distributed processing. Furthermore, we are often overwhelmed 
by a deluge of data assigned to one single frame system, which becomes simply too 
large to be handled numerically. In these cases it would be highly beneficial to split a 
large frame system into a set of (overlapping) much smaller systems, and being able 
to process locally within each sub-system effectively. 

A distributed frame theory relating to a set of local frame systems is clearly in 
demand. In this paper we develop a suitable theory based on fusion frames, which 
provides exactly the framework not only to model these applications but also to 
provide efficient algorithms with sufficient robustness. 

1.1. Applications under Distributed Processing Requirements. A vari- 
ety of applications require distributed processing. Among them there are, for instance, 
wireless sensor networks ,2.Q. . geophones in geophysics measurements and studies |14j . 
the physiological structure of ear and hearing systems |53]. To understand the na- 
ture, the constraints, and related problems of these applications, let us elaborate a 
bit further on the example of wireless sensor networks. 
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In wireless sensor networks, sensors of limited capacity and power are spread in 
an area sometimes as large as an entire forest to measure the temperature, sound, 
vibration, pressure, motion and/or pollutants. In some applications, wireless sensors 
are placed in a geographical area to detect and characterize chemical, biological, 
radiological, and nuclear material. Such a sensor system is typically redundant, and 
there is no orthogonality among sensors, therefore each sensor functions as a frame 
element in the system. Due to practical and cost reasons, most sensors employed in 
such applications have severe constraints in their processing power and transmission 
bandwidth. They often have strictly metered power supply as well. Consequently, 
a typical large sensor network necessarily divides the network into redundant sub- 
networks - forming a set of subspaces. The primary goal is to have local measurements 
transmitted to a local sub-station within a subspace for a subspace combining. An 
entire sensor system in such applications could have a number of such local processing 
centers. They function as relay stations, and have the gathered information further 
submitted to a central processing station for final assembly. 

In such applications, distributed/local processing is built in the problem formula- 
tion. A staged processing structure is prescribed. We will have to be able to process 
the information stage by stage from local information and to eventually fuse them 
together at the central station. We see therefore that a mechanism of coherently 
collecting sub-station/subspace information is required. 

Meantime, due to the often-time unpredictable nature of geographical factors, 
certain local sensor systems are less reliable than others. While facing the task of 
combining local subspace information coherently, one has also reasons to consider 
weighting the more reliable sets of substation information more than suspected less 
reliable ones. Consequently, the coherent combination mechanism we just saw as 
necessary often requires a weighted structure as well. We will show that fusion frame 
systems are created to fit such weighted and coherent fusion needs. 

1.2. Parallel Processing of Large Frame Systems. In case that a frame 
system is simply too large to handle effectively from the numerical stand point, there 
are needs to divide the large system into smaller and parallel ones. Like many parallel 
processing mechanisms, one may consider splitting the large system into multiple small 
systems for simpler and parallel processing. Evidently, the subdivision mechanism 
must take into consideration a coherent combination after the subsystem processing. 
To make the subdivision mechanism more robust, one may not want to (sometimes it 
is also impossible to) split the large system in an independent or orthogonal fashion. 
Such a splitting and then a coherent combination must produce precisely the original 
result if the system were to be processed globally. 

Fusion frame systems are created to fit such needs as well. Weighted coherent 
combination of subsystems (as provided by fusion frame theory) is also useful in such 
applications where losses of some subsystem information occur. Sometimes, weighted 
coherent combination is also useful from an efficient and approximation point of view. 
Some approaches such as the domain decomposition method '29? also use coherent 
combinations. However the fusion frame theory will provide a much more flexible 
framework that also takes local frames into account. 

1.3. Fusion frames. In this article, we are interested in weighted sequences of 
subspaces with controlled "overlaps" . Each subspace is equipped with local frames 
aiming at the development of a framework for the applications discussed above. In 
[TT] , two of the authors studied redundant subspaces for the purpose of easing the con- 
struction of frames by building them locally in (redundant) subspaces and then piec- 
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ing the local frames together by employing a special structure of the set of subspaces. 
This was referred to as a frame of subspaces. Related approaches were undertaken by 
Aldroubi, Cabrelli, and Molter and Fornasier A similar idea was also used 

by Aldroubi and Groching in a quite different context in Moreover, Bodmann, 
Kribs, and Paulsen employed Parseval frames of subspaces for optimal transmission 
of quantum states 0- Some further results on the theory from m can be found in 
[3, and an extension was derived by Sun [271 128j . however without any possibility of 
equipping the subspaces with an underlying structure. 

We will employ some parts of this theory and set it into the context studied in 
this paper. The structure of the overlapping subspaces, e.g., which relate to the sub- 
networks of a wireless sensor network, will be modeled by employing the notion from 
To avoid confusion with the long existing term "frame for subspaces" and to 
emphasize the fact that this mathematical object will provide a framework to fuse 
data in the subspaces, we decided to coin it fusion frame in this context. In our 
situation it will become also essential to view a fusion frame together with a set of 
local frames for its subspaces, in which case we will speak of a fusion frame system. 

We observe that fusion frames contain conventional frames as a special case. This 
theory goes thereby "beyond frame theory" . It turns out that the fusion frame theory 
is in fact much more delicate due to complicated relations between the structure of 
the sequence of weighted subspaces and the local frames in the subspaces and due to 
the extreme sensitivity with respect to changes of the weights. 

Our main motivation is to study fusion frame systems with respect to their recon- 
struction properties to not only provide a comprehensive model for applications which 
require distributed processing and which employ a distributed structure due to com- 
plexity reasons, but also to build efficient algorithms for fusion and reconstruction. 
We provide a general reconstruction formula by employing a so-called fusion frame 
operator, derive a variety of ways to fuse/reconstruct depending on the ability of the 
application to process off-line or only in real time, and present an iterative algorithm. 
Since we are also concerned with applications having the choice between distributed 
and centralized reconstruction we further show that in very special cases those re- 
constructions are in fact performed by employing the same set of vectors, thereby 
presenting situations where distributed reconstruction demonstrates the same behav- 
ior as centralized reconstruction, e.g., with respect to noise. 

As discussed above, sensor networks in particular suffer significantly from distur- 
bances of individual sensors or even whole sub- networks in the form of, e.g., natural 
forces. This led us to the study of stability of fusion frame systems not only under 
perturbations of the subspaces themselves, but even more of the local frame vectors. 
In order to describe the properties of the affected sensor network explicitly, we present 
several results which, in particular, give precise estimates for the changes of certain 
properties of fusion frame systems. 

1.4. Contents. The organization of this article is as follows. In Sectional the 
definition of fusion frames and fusion frame systems and their fundamental charac- 
terization will be given. Examples of fusion frames are presented, and connections 
of fusion frames with conventional frames will be discussed. In Section ^ several 
fusion frame reconstructions are presented. These are the coherent combinations we 
discussed earlier. Both operator theoretical and its matrix representation are consid- 
ered. An iterative fusion reconstruction is also constructed in this section. Section 
iniis devoted to the robustness of fusion frames, in which the analysis of stability of 
fusion frame systems to perturbations is extensively carried out. Conclusion remarks 
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and application discussions are the subjects of the last section. 

2. Review of Frames and Notation. A sequence T ~ in a Hilbert 
space 7i is a frame for 7i, if there exist < A < B < oo {lower and upper frame 
bounds) such that 

A\\fr<J2\{fJ.)\'<B\\ff for all /GH. (2.1) 

The representation space associated with a frame is (■2{I)- In order to analyze a signal 
f € Ti., i.e., to map it into the representation space, the analysis operator Tyr ; Ji — + 
£2(1) given by Tjrf = {(/, fi)}i^i is applied. The associated synthesis operator, which 
provides a mapping from the representation space to H, is defined to be the adjoint 
operator : £2(1) 'H which can be computed to be T^{{ci\iizj) = X^ig/ Ci/i- By 
composing Tjr and we obtain the frame operator 

Sr-.H^n, Srf = T*T^f = Y,{fj^)f^. 

Whenever = {/ijig/ is a frame, we know that there exists at least one dual 
frame {fi}iei satisfying 

/ = E</'/»)/^ = E</'^')/' for all /G?^. (2.2) 

When is a redundant (inexact) frame, there exist infinitely many dual frames {fi}i£i 
- which can even be characterized |21| - in which the canonical dual frame defined by 
{S^^fi}i(=i is the one having the least square property among all dual frames {fi}iei- 
That is, for all feH, E.e/ < E.e/ 

Of particular interest are A-tight frames, i.e., if the frame bounds can be chosen 
as A = _B in the frame definition H2.1|) . Provided H2.1|l holds with A ~ B — 1, we 
call a Parseval frame. The advantage of working with these frames can be clearly 
seen by considering the reconstruction formula H2.2II . In these cases the canonical 
dual frame equals {-jfi}i£i, and hence we obtain / = ^T^Tjrf for each f G H, i.e., 
we can employ the frame elements for both the analysis and the synthesis. There 
exist many procedures to construct tight or Parseval frames (cf. O^]). However, 
Parseval frames with special properties are usually particularly difficult to construct, 
see, e.g., 

For more details about the theory and applications of frames we refer the reader 
to the books by Christensen Daubechies and Mallat 22 . 

3. Fusion Frames. In this section, the notion of a fusion frame and a fusion 
frame system is introduced. Discussions of the notion of redundancy for fusion frames 
is also provided. We will put our focus on the structure of the fusion frame operator 
and its connection with the fusion frame bounds, the reason being that the fusion 
frame operator will become essential for studying distributed fusion/reconstruction in 
Sectional Finally, the fact that our theory goes "beyond frame theory" in the sense 
that conventional frames are a special case of fusion frames will be discussed, and 
how much more sophisticated the theory of fusion frames turns out to be will also be 
highlighted. 
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3.1. Definition and Basic Properties. We will start by stating the definition 
of a fusion frame. 

Definition 3.1. Let I be some index set, let {M^ijie/ ^6 family of closed 
subspaces in H., and let {vijig/ be a family of weights, i.e., Vi > for all i Cz I . Then 
{(Wi, Vi)}i^i is a fusion frame, if there exist constants < C < D < oo such that 

qi/f < E II '^M^. (/)!!' ^ ^11/11' / e (3.1) 

iei 

where irwi is the orthogonal projection onto the subspace Wi. We call C and D the 
fusion frame bounds. The family {(Wi,Wi)}ig/ is called a C-tight fusion frame, if 
in \3. 1\) the constants C and D can be chosen so that C — D, a Parseval fusion 
frame provided that C = D = 1 and an ortlionormal fusion basis if Ti. — ©jg/ Wi. 
If {{Wi,Vi)}i^i possesses an upper fusion frame bound, but not necessarily a lower 
bound, we call it a Bessel fusion sequence with Bessel fusion bound D. 

Often it will become essential to consider a fusion frame together with a set of 
local frames for its subspaces. In this case we will speak of a fusion frame system. 

Definition 3.2. Let {{Wi,v.i)\ii^i be a fusion frame for Ti, and let {fij}j^j. i^j 
be a frame for Wi for each i <E I . Then we call {{Wi,Vi, {/ijjjgjjjie/ a fusion frame 
system for Ti.. C and D are the associated fusion frame bounds, if they are the fusion 
frame bounds for {{Wi, Wi)}ie/; and A and B are the local frame bounds, if these are 
the common frame bounds for the local frames {fij}jeJi for each i d I. A collection 
of dual frames {fij}jeJi, i ^ I associated with the local frames will be called local dual 
frames. 

For a fusion frame system we have the following intriguing relation between prop- 
erties of the associated fusion frame and the sequence consisting of all local frame 
vectors, in this sense it provides a link between local and global properties. For the 
proof we refer to ^2 Thm. 3.2]. 

Theorem 3.3. For each i € I, let Vi > 0, let Wi be a closed subspace ofTL, and 
1st {fij}jeJi a frame for Wi with frame bounds Ai and Bi. Suppose that 

< A = inf < sup Bi ^ B < oo. 

iei 

Then the following conditions are equivalent. 

(i) {(Wi, Vi)}iei is a fusion frame for Ti. 

(ii) {vifij}je.h,tei is a frame for H. 

In particular, if {(Wi,Vi, {fij}j,= j-)}i,=i is a fusion frame system for Ti. with fusion 
frame bounds C and D, then {vifij}jeJi.i^i is a frame for Ti with frame bounds AC 
and BD. Also if {vifij}i^ij^j^ is a frame for Ti with frame bounds C and D, then 
{(Wi, Ui, {/y}jgj.)}ig/ is a fusion frame system for Ti with fusion frame bounds ^ 
and ^. 

Tight frames play a vital role in frame theory due to the fact that they provide easy 
reconstruction formulas, and also tight fusion frames will turn out to be particularly 
useful for distributed reconstruction (cf. Section^. The previous theorem implies 
the following relation between tight fusion frames and tightness of the collection of 
the local frames in a fusion frame system. 

Corollary 3.4. For each i d I, let Vi > 0, let Wi be a closed subspace of Ti, 
and let {fij}jeJi be a Parseval frame for Wi. Further, let C be a constant. Then the 
following conditions are equivalent. 

(i) {(Wi, Vi)}i^i is a C-tight fusion frame for Ti. 
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(ii) {vifij}je.h,i<£i «s a C -tight frame for H. 

By employing this result, the redundancy of a finite C-tight fusion frame can be 
made precise in terms of the fusion frame bound. 

Proposition 3.5. Let {(Wi,Vi)}^^i be a C-tight fusion frame for H. Then we 
have 



C 



dimTi 



Proof. Let {ey }^™^' be an orthonormal basis for Wi for each 1 < i < n. By 
Corollarv l3.4l the sequence {wieij}"^;^ j™^' is a C-tight frame for Ti. Employing ^1 
Sec. 2.3] yields that 

dim Ti dim Ti 

□ 

In this sense, we can interpret the frame bound C as the redundancy of the tight 
fusion frame {(Wi, Wi)}"^^. 

To enlighten the definitions let us consider the following example. 

Example 3.6. Since almost all applications require a finite model for their 
numerical treatment, we restrict ourselfs to a finite dimensional space in this example. 
Suppose {fn}n=i ^ frame for M.^ with frame bounds A, B. Now we split {1, . . . , N} 
into K sets Ji, . . . , Jk, and define Wi = span{/„}„gj., 1 < i < K . Since in the finite- 
dimensional situation each finite set of vectors forms a frame, in particular {fn}ne.h 
is a frame for Wi for each \ < i < K . Let C and Z? be a common lower and upper 
frame bound, respectively. Theorem 13.31 now impfies that {(Wi, 1, {/njne Ji)}^i is a 
fusion frame system with fusion frame bounds §, -j- Suppose that by weighting the 
subspaces we can make {(Wi, ■yi, {/„}„gj.)}^j a tight frame, then the fusion frame 
{{W,,v{)}ti has redundancy (^fli vUJi)/^. 

3.2. Fusion Frame Operator. In frame theory an input signal is represented 
by a collection of scalar coefficients that measure the projection of that signal onto 
each frame vector. The representation space employed in this theory equals 
However, in fusion frame theory an input signal is represented by a collection of 
vector coefficients that represent the projection (not just the projection energy) onto 
each subspace. Therefore the representation space employed in this setting is 

( $^©W, J = {{f^}^eI\f^ G W, and {II/.IIW e f{I)}. 
Kiel / £2 

Let W = {(Wi,Vi)}ig/ be a fusion frame for H. In order to map a signal to the 
representation space, i.e., to analyze it, the analysis operator Tvv is employed, which 
is defined by 



w ■ n-^ I ^©Wj j with Tw(/) = {v^'^w^{f)}t<£I■ 
\iei ) 

It can easily be shown that the synthesis operator T-^, which is defined to be the 
adjoint operator, is given by 



^ ®W, j -^n with Twif) - J2 / = e [ ^ ®W, j 

J e„ iei Kiel / , 
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Now we can give the definition of a fusion frame operator. The fusion frame operator 
Sw for W = {{Wi, Vi)}i(zj is defined by 

Interestingly, a fusion frame operator shows results similar to a frame operator 
concerning invertibility. For the proof of the following result we refer to [111 Prop. 
3.16]. 

Proposition 3.7. Let {(Wi,Vi)}i£i be a fusion frame for TL with fusion frame 
bounds C and D . Then the associated fusion frame operator S'w is a positive and 
invertible operator on Ti. with C Id < S'w < D Id. 

3.2.1. Fusion frame operator in terms of local frames. For the purpose 
of distributed fusion/reconstruction (see Section^, employing a fusion frame system 
the fusion frame operator will indeed become essential. More precisely, the inverse of 
the fusion frame operator will be employed. Therefore, a further investigation of the 
fusion frame operator computationally is helpful. 

We observe that the fusion frame operator can be expressed in terms of local 
frame operators as follows: 

Proposition 3.8. Let {{Wi,Vi,J-'i — {fij}jeJi)}iei be a fusion frame system for 
Ti., and let Ti — {fij}jeJi, i € I be associated local dual frames. Then the associated 
fusion frame operator Sw can be written as 

Proof For all f eli., 

swf = E ^ E E (/' /^^■)^- = E E (/' 

Applying the definition of the analysis operators Tjr. , and the associated synthesis 
operators (see Section |21l, the result follows immediately from here. □ 

3.2.2. Matrix representation of the fusion frame operator. For computa- 
tional needs, let us further consider the fusion frame operator in finite frame settings, 
where the fusion frame operator will become the sum of (weighted) matrices of each 
subspace frame operator (see also Example 13. 6|l . 

Let Fi be the frame matrices formed by frame vectors {/ijjjgji in the column- 
by-column format 

Pi = {fil fi2 ■ ■ ■ fiji) ■ 

Similarly, let Fi be defined in the same way by the dual frame {fij}j^j.. Then the 
fusion frame operator associated with finite frames has the expression 

5W = E = E 

i£l i£l 

where stands for the Hermitian transpose of a matrix M . Therefore, the evalua- 
tion of the fusion frame operator and the inverse fusion frame operator in finite frame 
settings are quite straightforward. In practical applications, this will turn out to be 
very convenient. 
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3.3. Fusion Frame Bounds. Our first aim is to establisli a connection between 
the fusion frame bounds and the norm of the fusion frame operator. We achieve this 
by first showing that the boundedness of the associated fusion frame operator is 
equivalent to the weighted sequence of closed subspaces satisfying the fusion frame 
property. 

Proposition 3.9. Let {Wi}i^i be closed subspaces in Ti., let {vi\i^i be positive 
numbers, and let Syy denote the fusion frame operator associated with {{Wi, 

(i) If Syv < D Id, then {{Wi,Vi)}i^j is a Bessel fusion sequence with bound D. 

(ii) If Sw > C Id, then {{Wi,Vi)}i^i possesses the lower fusion frame bound C. 
Proof, (i). Let Tyy denote the analysis operator associated with {{Wi,Vi)}i^i. 

Since Sw = Tw'^w ^^'^ hence ||Tw|P = H^wH, for any f gH we obtain 

E^'li'^^. = WT^^-fW ^ WTwW'Wfr = WSwWWfW < DWfW^ 

iei 

(ii). For all f G H, we have 

rw/ip - {T^Ty,f,f) = (Swf, f) = {sif,sif) = wsiff > c\\f\\\ 

□ 

Theorem 3.10. Let {Wi}i£i be closed subspaces in H, let {wijig/ be positive 
numbers, and let Sw denote the fusion frame operator associated with {{Wi,Vi)}i£i . 
Then the following conditions are equivalent. 

(i) {(Wi,Wi)}ig/ is a fusion frame with fusion frame bounds C and D. 

(ii) We have C Id < Sw < D Id. 

Moreover, the fusion frame bounds are ||S'wll 0,''^^ ||S'yy'^||. 
Proof, (i) => (ii). This is implied by Proposition 13. 71 
(ii) ^ (i). This follows from Proposition |^| □ 

Next we will study the behavior of the fusion frame operator under applying a 
self-adjoint and invertible operator to the set of subspaces. In Subsection 13.41 this 
result will reveal essential differences between fusion frame theory and classical frame 
theory. 

Proposition 3.11. Let {{Wi,Vi)\i^i be a fusion frame for Ti. with associated 
fusion frame operator Sw, d'^d let T be a self-adjoint and invertible operator on TL. 
Then {(TW^, i'i)}ig/ is a fusion frame for H with fusion frame operator TSwT~^. 

In particular, {{S^Wi,Vi)}ifzi anrf {(5'vy W^, Wi)}.;^/ both possess the fusion frame 
operator Swj ^.i^d, hence are fusion frames with the same fusion frame bounds as W. 

Proof. We recall the following basic fact from frame theory: Provided that = 
{fi}iei forms a frame for a closed subspace W oi H with frame operator Sj^, then 
{Tfi}i£i is a frame for TW with frame operator TSj^T, since 

Y^if, Tf,)Tf, = T [Y^^Tf, /,)/,) - TS:rTf. 
iei Kiel J 

Hence the dual frame of {r/i}ig/ is {T^^Syr^ fi}i,=i. 

For each i G I, let Ti — {fij}jeJi be a frame for Wi. Then for all / G we 
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compute 



iei iei \jeJi 



) 



Vie/ 



= TSwT-^f. 



The second part follows from here immediately by employing Theorem 13. 101 □ 

This result shows in particular that the associated fusion frame is invariant under 
the application of a self-adjoint and invertible operator, which commutes with the 
fusion frame operator, to the set of subspaces. 

3.4. Beyond Frame Theory. Interestingly, frames can be shown to be a special 
case of fusion frames in a particular sense with the natural meaning of the fusion frame 
bounds and the fusion frame operator. We will make this precise in the following 
proposition. 

Proposition 3.12. Let T — he a frame for H with frame bounds A, B. 

Then {(span{/i}, is a fusion frame for Ti. with fusion frame bounds A, B 

and fusion frame operator Syr . 

Proof. Observe that for any f E H, 



Hence the fusion frame operator for {(span{/i}, equals Sjr. Since 

possesses the frame bounds A and B, it follows that Aid < Sjr < Bid. Now Theorem 
13.101 implies that {(span{/i}, is a fusion frame for Ti. with fusion frame 

bounds A and B. □ 

This result seems to indicate that fusion frame theory is "just" a generalization 
of frame theory. However, in the following remark we will enlighten the much more 
delicate behavior of fusion frames. A variety of further essential differences will be 
revealed by the results in the following sections. 

Remark 3.13. To demonstrate the much more rich behavior of fusion frames in 
contrast to frames, we consider a frame J- = {fi}i£i for Ti. with frame bounds A, B. 
Proposition 13 . 1 21 implies that 



is a fusion frame for Ti with fusion frame bounds A, B and fusion frame operator Syr. 
Since {S^^ fi}i^i is the canonical dual frame for {fi}iei with bounds B^^,A^^ and 
frame operator S^^, also 



is a fusion frame, but now with fusion frame bounds B~^,A~^ and fusion frame 
operator S^^. Now it is possible to change the associated weights to "move" the 
associated fusion frame operator back into the range of Aid and Bid. This is done 
by applying Proposition 13 . 1 II to {(span{/i}, which yields that 



E ll/»ll''^span{/,}(/) = E m^m = E</' ^ ^^f- 




{(span{/,;},|l/,|l)},g, 



{(span{5^V»},II^^V.||)he/ 



{(5^ispan{/,}, = {(span{5^ 
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is a fusion frame with bounds A, B and fusion frame operator Syr. This observation 
not only reveals how much more sensitive fusion frames behave, but also indicates 
how critical the selection of the weights can be. 

Our observation is based on the fact that under the application of a self-adjoint 
and invertible operator T to both a frame and a fusion frame, the frame operator 
changes to TSj^T, however the fusion frame operator changes to TSwT~^. Note that 
this fact for frames ensures that applying Sjr^ to the frame vectors yields a frame with 
frame operator S^^. However, the analog formula for fusion frames seems to make 
it almost impossible to construct a new fusion frame having as a fusion frame 
operator. We like therefore to state this as an open question. 

4. Distributed Fusion/Reconstruction. Given a large set of data, some ap- 
plications such as certain data fusion problems jSO] require processing the data first 
locally by employing a frame structure, and then fusing the (computed) subspace in- 
formation globally. This procedure is called distributed fusion, and obviously, the sec- 
ond step can be modeled by employing the framework of fusion frame systems. If the 
initial data comes from a decomposition of a signal with respect to a global frame such 
as in sensor networks problems |2(Jj . and the task consists in precisely reconstructing 
the initial signal via the procedure mentioned above, we speak of distributed recon- 
struction. In this case we sometimes do have the choice of whether either performing 
distributed or centralized reconstruction, an issue that will be further elaborated in 
this section. 

We shall first analyze the different distributed fusion procedures depending on 
whether it is necessary to perform real time operations or whether it is possible to 
compute certain operations off-line. We will then present an iterative algorithm for the 
computations in these procedures. Finally, we discuss several aspects of distributed 
reconstruction versus centralized reconstruction, in particular concerning the set of 
vectors employed to perform the reconstruction. 

4.1. Distributed Fusion Processing. The first fundamental observation we 
make consists of the fact that distributed fusion processing is feasible in an elegant 
way by employing the inverse fusion frame operator. 

Proposition 4.1. Let {{Wi,Vi)}i^i be a fusion frame for Ti. with fusion frame 
operator Syy and fusion frame bounds C and D . Then we have the reconstruction 
formula 

f^Y. ^^Sy,'7^w, if) for all f e H. 

Proof. By Proposition 13. 71 for all f E H we have 

/ = SyJSwf = ^ v^SyJ-^^w, if)- 
iei 

□ 

The fusion frame theory in fact provides two different approaches for distributed 
fusion procedures. For this, let {{Wi,Vi, {fij}j^j.)}i^i be a fusion frame system for 
TC, and let {fij}jeji, i £ I he associated local dual frames. 

One distributed fusion procedure is from the local projections of each subspace: 

/ - E ^'S^'^wJ = E ''^^w I E (/' /»^)/*^ I fo-^ all / ^ ^- (4.1) 

iei iei \jeJi I 
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In this procedure, the local reconstruction takes place first in each subspace Wi, and 
the inverse fusion frame is applied to each local reconstruction and combined together. 

Another form of distributed fusion actually acts like a global reconstruction if the 
coefficients of signal/function decompositions are available: 

/ = E E (/' /'^) (Swk) for aU fen. (4.2) 

iei jeJi 

The difference in this fusion procedure compared with global frame reconstruction 
lies in the fact that the (global) dual frame {S^fij} is first calculated at the local 
level, and then fused into the global dual frame by applying the inverse fusion frame 
operator. This makes the evaluation of (global) duals much more efficient. 

Remark 4.2. Depending on applications, some may require the fusion procedure 
via H4.1|l such as in sensor networks 20 , and geophones in geophysics measurements 
|14|. whereas some may allow for fusion process via (|4.2() such as parallel processing 
of large frame systems. Let us examine the orders of computation of the fusion pro- 
cedures H4.1|l and 14.2|l . respectively. Besides the operation of in both equations, 
both fusion procedures have the same number of multiplications. However, (14.1(1 typ- 
ically has less (but real time) inverse fusion frame operations. Specifically, (|4.1|l has 
|/| operations of over |/| local reconstructions. On the other hand, 14.211 requires 
Sig/I'^il operations of over local dual frames {/^ j^gj^^ig/, which is typically 
much larger than the |/| operations in 1(4. l|l . It is nevertheless equally important to 
point out that the much larger operation requirement in ((4.2(1 can be carried out 
"off-line", which often-times can be advantageous. 

4.2. Iterative Reconstruction. Fusion frame reconstruction can be carried 
out iteratively as well, just like in frame reconstructions The specific mechanisms 
can also be divided in two different ways, depending on whether a local reconstruction 
actually takes place or not as given in ((4.1() or ((4.2() . 

The first way we present refers to the distributed fusion procedure given by ((4.1(1 . 

Proposition 4.3. Let {{Wi,Vi)}i^i be a fusion frame in Ti, with fusion frame 
operator and fusion frame hounds C , D. Further, let f G Tl, and define the 
sequence (/„)„gNo by 

fo, n = 0, 

\fn -1 + cTS'^w(/ - /n-i), n > 1. 

Then we have f — lim„^oo fn with the error estimate 



Proof. Employing the construction of the sequence (/rOneNoi for each n e N we 
obtain 

f - fn ^ f - fn-l - (g -\- D ^^^'^ ~ fn-l) = \J ^ C + D ^^J ^ fn~l)- 



Interating this argument yields 

-Jr- 

C + D 



f-.fn = [l-TT^Sw] (/-/o). (4.3) 
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Now we use Proposition 13 . 71 to conclude that 
In a similar way, we can show that 



Combining these two estimates yields 



2 „ „ D-C ^ 
SwW < TTTT^ < 1. 



C + D - D + C 
Applying this estimate to (|4.3|l . we finally obtain 

2 / D — C\ " 

||/-/„||2<||/-^^^W||"||/-/0|| < [-^:^) 11/11- 

□ 

Thus every f G Ti, can be reconstructed from the fusion frame coefficients Tw (/) = 
{viTTWi (/)}ie/) since Swf does only require the knowledge of those and of the sequence 
of weights {vi}i£i. In each iteration, Swf is always known. The only significant 
computation is Swfn, which can be easily carried out using the local frame structure 
of each subspace. 

We remark that an application of the Chebyshev method or the conjugate gradi- 
ent method as done by Grochenig 18^ for the frame algorithm should lead to faster 
convergence. 

Finally, we discuss an interactive way to compute the distributed fusion procedure 
given by l|4.2|l . 

Remark 4.4. In some applications, if the local measurements/local frame coef- 
ficients are preserved, the final reconstruction can also be done in a "global" fashion 
with distributed evaluation of (global) duals through local dual frames {fij}je,Ji ,iGl 
associated with a fusion frame system {{Wi,Vi, {fij}j(EJi)}iei- 

This follows from equation 14.2|l . in which the evaluation of local duals {fij}jeJi 
is carried out distributively, and the evaluation of {S^ fij}j^j-^i^i can be carried out 
iteratively as described in Proposition 14.31 

4.3. Distributed Reconstruction and (Global) Dual Frames. The pur- 
pose of this section is to study the sequence of vectors employed for distributed re- 
construction and to compare distributed with centralized reconstruction. For this, let 
{(Wi, Vi, {fij}ji£ji)}iei be a fusion frame system for Ti. with local frame bounds A, B, 
and let {/ijjjeJij i € / be associated local dual frames. Since, by Theorem 13.31 the 
sequence T = {vifij}j^j.^i^i is a frame for 7i, we might consider the situation that 
we are given the (global) frame coefficients {(/, Vifij)}j^j.^i^i of a signal f £H. For 
some applications, which do not enforce distributed reconstruction, we might have 
two ways to reconstruct /. The (global) dual frame {S^^Vifij}j£j-^i^i could be used 
to perform centralized reconstruction^ i.e., to compute 

/ XI XI ^»/y ) (S^^'^^h) ■ 
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Or, in order to reduce the complexity, we might employ the associated fusion frame 
operator S'yy to perform distributed reconstruction, and obtain (compare H4.2|l ') 



In the sequel we will discuss the difference between the sequences {S^^Vifij}j^j.^i^j 
and {S^Vifijjjej^^iei in more detail. 

Our first result shows that indeed {S^Vifij}j(:j.^i£i is a dual frame for !F, but 
not necessarily the canonical dual frame. 

Proposition 4.5. Let {{Wi,Vi,{fij}j^j.)}i(=i be a fusion frame system for H 
with associated fusion frame operator S'w, local frame bounds A, B and local dual 
frames {fij}j<£Ji, i € I. Then {S^Vifij}j(zj.j^i is a dual frame for the frame 
{'Vifij}jeJi,i<£i ■ 

Proof. First we note that {vifij}j^j.^i^j indeed forms a frame by Theorem 13.31 
Employing the self-adjointness of Sw fProDOsition l3.7|l . we have for all / G 

iei jeJi iei jeJi 

iei jeJi 
iei 

= SwiSyJ'f) 

which finishes the proof. □ 

It is interesting to observe that a "dual" relation also holds. We wish to mention 
that this property does not have quite the same correspondence in conventional frames 
as well. 

Proposition 4.6. Let {iWi,Vi,{fij}je.ii)}iei be a fusion frame system for H 
with associated fusion frame operator Sw, local frame bounds A, B and local dual 
frames {fijjjej,, i € /. Then {wj/yl^g/jg is a frame for H and {S^^^Vifi-jji^ijeJ. 
is a dual frame for it. 

Proof. The fact that {I'i/yjie/.je is a frame for Ti. follows again from Theorem 
13.31 since each local frame {fij}jeJn i & I has frame bounds A~^, B"^ and thus 
possess a common lower and upper bound. 

Now using the fact that {vifij}i^j,j^j. and {S^Vifij}i^jj^j. are a pair of dual 
frames of H by Proposition 14.51 we have for all / G Ti, 

f ^^^{f' Vif,j)S^vJ,j 
iei jeJi 

\iei jeJi 

\iei je.ii 
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U 

In order to compare distributed reconstruction with centralized reconstruction, it 
is essential to understand when {S^Vifij}j(zj.j^j equals the canonical dual frame 
of the frame = {vifij}j(zj.j^i (compare also (|4.1|l '). since in these particular cases 
distributed and centralized reconstruction coincide. In general, this certainly need 
not be the case due to the observation that if we have a Parseval fusion frame, then 

with the Sjr. being the local frame operators, and hence due to the occurring cross 
terms, 

However, the following results give some special cases in which distributed and cen- 
tralized reconstruction indeed coincide. 

Proposition 4.7. Let {{Wi,Vi,J^i — {fij}jeJi)}iei be a fusion frame system for 
TL with associated fusion frame operator Sy\; , local frame bounds A, B and local dual 
frames {fij}ji£ji, i € I. If {{Wi,Vi)}i,£i is an orthogonal fusion basis or {fij}jeJi 
a Parseval frame sequence for all i Cz I , then {S^ViS^^ fij}j^j._i£i is the canonical 
dual frame of the frame {vifij}ji£j-_i£i. 

Proof. By Theorem 13.31 the sequence J- = {vifij}j£j-^i,=i forms a frame for Ti., 
and we denote its frame operator by Sj^. If {(W;, fi)}ig7 is an orthogonal fusion basis, 
then Sy\! = Id and Syr = X^ig/ , and hence S^^ = X^ig/ (BS^^TTWi ■ Provided 

that {fij}jeJi is a Parseval frame sequence for all i € I, we have Sj^- = Id for all 
i G I, and we further obtain for all / e H, 

Swf = '"f'^w, (/) = X! X! = ^^Ff- 

iei iei jeJi 

In both cases the claim follows immediately from here. □ 

Finally, we would like to point out a surprising fact concerning the situation of 
having Parseval frames spanning the subspaces of a fusion frame, which arises from 
this result. 

Remark 4.8. Let {{Wi,Vi,!Fi = {/y IjejJIie/ be a fusion frame system for H 
with associated fusion frame operator S'yy, and let J^i be Parseval frames for all i £ I . 
By the previous result, the operator Sw is independent of the choice of the Parseval 
frame, since Sw always equals the frame operator of the frame ig/. The 

intuitive reason for this is that provided we take Parseval frames for the subspaces, 
the frame property of the total collection of frame elements completely mirrors the 
behavior of the fusion frame. 

Furthermore, we would like to briefly mention the impact of Proposition 14.71 on 
the noise reduction under distributed and centralized reconstruction of this result. 

Remark 4.9. In [23 Rozell and Johnson studied noise reduction under dis- 
tributed reconstruction versus centralized reconstruction in finite dimensional Hilbcrt 
spaces. They used additive zero mean white noise, which they added to all frame 
coefficients with respect to the frame T = {vifij}j^j-j,=i and derived bounds for the 
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mean square error for distributed and centralized reconstruction. By numerical simu- 
lation they showed that randomly adding elements to a fusion frame not only improves 
the mean square error for distributed reconstruction, but even forces it to converge to 
the mean square error for centralized reconstruction. Goyal, Vetterli, and Thao |17| 
proved that a frame becomes asymptotically tight if elements are randomly added. 
By employing this result it was argued in |2Hj that some very restrictive conditions 
on the local frames and the fusion fame structure might be fulfilled in the limit. 

Applying Proposition 14.71 now gives a broader picture for this intriguing phe- 
nomenon. By |17| . the local frames become asymptotically tight with the same frame 
bound, say A, as more and more random elements are added. Hence, by Proposition 
14.71 the sequence {S^ViS^^^^fij}j^j-jizi, which is employed for the distributed re- 
construction as outlined in (j4.1|l (except the constant ) , converges to the canonical 
dual frame of the frame {vifij}j^j-^i^i, which is used for the centralized reconstruc- 
tion. Taking the fact into account that the additional multiplicative constant does 
not play any role concerning the noise reduction ability, reveals a reason for the phe- 
nomenon described in [23j and does not require restrictive conditions. 

5. Robustness of Fusion Frame Systems. In this section we analyze the 
stability of fusion frame systems under perturbations of both the subspaces which 
constitute a fusion frame and the local frame vectors contained in the subspaces. The 
reason for this is that on the one hand, for instance, several complete groups of geo- 
phones might be moved to a slightly different location to adjust for transmission 
conditions, and on the other hand, for instance, in wireless sensor networks the loca- 
tion of single sensors might be changed slightly due to the impact of natural forces 
pn|. Therefore it is essential to study the robustness of fusion frame systems under 
these two different impacts. 

Thus, with these practical aspects in mind, we proceed by first examining pertur- 
bations of the subspaces in Section 15.11 and secondly studying robustness of a fusion 
frame system under perturbations of the associated local frames in Section 

5.1. Perturbation of the Fusion Frame. First we would like to point out 
one fundamental problem with perturbations of fusion frames which is the cause of 
the serious technicalities in these results. Since the main ingredient in the definition 
of a fusion frame are the orthogonal projections onto a set of subspaces, it would 
be natural to consider perturbations of these projections. However, there is no such 
thing as a perturbation of a projection. This means, that if P and Q are projections 
on 7^, < Ai , A2 < 1 , and 

\\Pf - QfW < AillP/ll + A2IIQ/II for aU f e H, 

then it follows that P — Q. This can be easily seen by way of contradiction as follows: 
If P Q, then there exists a vector f G H so that / _L P{Ti), but also satisfying 
Qf 7^ (or vice- versa). This yields 

\\Pf - QfW = 110/11 < AiiiP/ii + A2II0/11 = A2||g/||, 

which is a contradiction. 

Therefore, we define (Ai, A2)-perturbations of sequences by employing the canon- 
ical Paley- Wiener-type definition: 

Definition 5.1. Let {Wi}ig/ and be closed subspaces in Ti, let 

be positive numbers, and Zei < Ai, A2 < 1 and e > 0. // 

Wi-KWi - Ti'^y.)/!! < Ai||7rvi/,/|| + A2||7r,;^,J|| +e||/|| for all f eU and i £ I, 
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then we say that {(Wj, ?;,)}jg/ is a (Ai, A2,£)-perturbation of {{Wi,Vi)}i^i. 

Employing this definition, we derive the following result about robustness of fusion 
frames under small perturbations of the associated subspaces. 

Proposition 5.2. Let {(Wi,Vi)}i^i be a fusion frame for H with bounds C, D. 
Choose < Ai < 1 and e > such that 

^l-X,)VC-e^^vij >0. 

Further, let {{Wi,Vi)}iei be a {Xi,X2,e) -perturbation of {{Wi,Vi)}iQi for some < 
A2 < 1. Then {{Wi,Vi)}i^i is a fusion frame with fusion frame bounds 



{l-X,)VC-e{j:,^jvf) 
I + A2 



1/2 



D(l + Ai)+£(E,e/^i 
1- A2 



Proof. We first prove the upper bound. For each f gH, we get 
1/2 



< 



< 



< 



1/2 



E ((1 + AiKlkm(/)ll + ^2V.h^^{f)\\+evm 
iei 

1 r 

iei 

(1+Ai) i^vUnwMW +A2 J2^n^wSm'- 



1/2 



1/2 



.iei 



liei 



1/2 . s 1/2 

1 1/2 . ^ 1/2 



Hence, 

(I-A2) 
which yields 



-,1/2 



12^ihwSf)\\' 

iiei 



< (1 + Ai) 



-1 1/2 



E^iii^^.(/)ii' 

Lie/ 



1/2 



Vie/ 



Lie/ 



1/2 



< 



^/^(l + Al)+£(Eie/^i) 
I-A2 



1/2 
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To prove the lower bound, for all / G we have 

-I 1/2 



> 



> 



1/2 



T.'''{\\''wm\-hwM)-rr^,im) 

lei 

E^' - Ai|kM..(/)|| - X2h^Sf)\\ - ^ll/ll) 

lei 

^ ((1 ~ AiK||7r;y.(/)|| - \2V^\^^m\ ~ evM\) 



1/2 



1/2 



> 



.lei 
This implies 



1/2 



lei 



1/2 



E^Vii/ii^ 



is/ 



1/2 



(1 + A2) 



lei 



1/2 



>(1-Ai) 



- 1/2 / \ 1/2 

E-'ii-^.(/)ii' -An-"] 11/11 

leJ J \iei / 

>{i-\,)Vc\\n-e(Y^vi^ 11/11 



which leads to 



E"'ii^H>.(/)ii 



iei 



1/2 



> 



(l-Ai)yC-e(E,6/«l 



,1/2- 



I + A2 



l/ll 



This completes the proof. □ 

We remark that a different perturbation result for fusion frames can be derived 
from [281 Thm. 3.1] by employing a different definition of perturbation. However, this 
would not lead to a result about robustness of fusion frame systems under disturbances 
of the local frames, which is what we are in particular aiming for. 

5.2. Perturbation of the Local Frames. The second fundamental problem 
with perturbing fusion frames locally is that a local perturbation cannot "see" the 
global structure of the fusion frame and therefore cannot adjust for it. 

For the notion of perturbations of sequences we employ the canonical Paley- 
Wiener-type definition (compare 0): 

Definition 5.3. Let and {fi\iei be sequences in H, and let < Ai, A2 < 

1. // 

II E«*(/* - ^ ^ill E°^/^ll + ^2|| E«^/^ll V all {aj,e/ e £^(1), 



iei 



iei 



iei 
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then we say that {fi}i^i is a (Ai, A2)-perturbation of {fi}i^i. 

First we derive properties of the relation between the two subspaces spanned by 
a sequence and its perturbed version. 

Proposition 5.4. Let {fi}i(=i be a frame sequence in TL, and let < Ai, A2 < 1. 
Suppose that {fi}i£i is a {Xi, X2) -perturbation of {fi}i,=i . Then is equivalent 

to In particular, we have 

\-y\\T.^^M ^ II E«»/'ll ^ ^^11 E«*/'ll f'^' e 

^ iei iei ^ iei 

and dim(spanjgj{/J) = dim(spanjgj{/i}). 

Set W = spanjgj{/i} and W = spauj^jj/i}. Then 



lkH^(7r^(/))|| > T— i - Ai-— ^ - A2 h^{f)\\ for all feH, 



1 - Ai _ ^ I + A2 
1 + A2 " Ai 

i.e., i/ Ai, A2 < ^, then tt^t is an isomorphism on Rng tTj^. 
Proof. The first part foUows from |13| . 
Let be the frame operator of = {fi}iei- For / £ H, 

II - /*)ii ^ ^iii J2^f^s^'f^)M + ^2ii E(/'^^'/^)/^ii 

iei iei iei 

^ iei iei 
Ai^ + A2) |k^(/)||. 

Employing this relation, we obtain 



= iiE</'^^'/')'^^(/»)ii 



iei 

> II E</' ^^'/^>'^1^(/0I1 - 11 E</' S~'f^hwif^ ~ m 

iei iei 

= II E(/' ^^'/^)/^ii - II E</' s^'f^)^w{fr - /Oil 

iei iei 



^:^llEa5^7.)/.ll - hw (^E(,/,^^7.)(/.-/,:)^ 



1 


- Ai 


1 


+ A2 


1 


-Ai 


1 


+ A2 


1 


-Ai 



> 



> — ^iiE(/'^^'^')^'ii - iiE</'^^'^')(/^ -/')ii 

iei iei 

>^IK^(/)II-(ai^ + A2) lk^(/)ll 



'"'^-Ai^-A2)||.^(/)||. 



I + A2 ^'^1-Ai ^"V"^ 
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It remains to observe that Ai, A2 < ^ implies — Ai ^^^ ^ — A 2 > 0. □ 

Remark 5.5. We wish to mention that in Proposition l5.4l we did not make use 
of the frame bounds of {fi}iei, but only of the constants Ai,A2 associated with the 
perturbation. Therefore it follows that our argument is symmetric in ttw and tt^^ and 
each permutation yields the same bounds. 

The following theorem gives a precise statement of how a perturbation of the 
local frames of a fusion frame system - which certainly results in a perturbation of 
the associated fusion frame - affects its fusion frame bounds, thereby in particular 
providing us with conditions under which the subspaces associated with perturbed 
local frames still constitute a fusion frame. Thus fusion frames systems are indeed 
robust not only against perturbations of the associated fusion frame fProposition l5.2|l , 
but even against perturbations of the local frames. 

Theorem 5.6. Let {{Wi,Vi, {fij}j^j.)}i^i be a fusion frame system for Ti. with 
fusion frame bounds C , D. Choose < Ai, A2 < 1 and e > such that 

For every i G /, let {fij}je.li be a {Xi, X2)- perturbation of {fij}j^j. and let Wi = 
spanj/ijljg J. . Then {{Wi, Vi)}i^i is a fusion frame for Ti. with fusion frame bounds 





2 




VC - e (Y.vf^'^^ 


and 











Proof. Fix i E I. Recalling Proposition l5.4l for all / e Ti we have 

hwm\'^hwMwM)\\' + \\ii'^wJ^wAf)f 

> (1 - y) lkm(/)f + \\{I-n^^)nwm\'- 

Hence, 

\\ii-^w.)^wAf)r<'-^hw.if)r- 

Employing Remark 15.51 Proposition 15 . 41 also yields 

Collecting the estimates derived above, for any f E H we obtain 
IKttw, - T^w,)if)f = ((ttw, - 7r^^)(/), (ttw, - 7r^J(/)) 

< - n^){i^wAf)) + {I - ^wM^^s^mwfw 

<||(/-7r^^)(^H'.(/))||||/|| + ||(/-7r^J(7r^^(/))||||/|| 
<ylkH.,(/)||||,/|| + ylk^,(/)||||/|| 

The theorem now follows from Proposition 15. 21 □ 
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6. Discussion. Fusion frames and fusion frame systems are natural extensions 

of the theory of frames. Wc have seen pressing needs of such notions and matlic- 
matical means in a variety of applications raging from sensor networks in geophysics, 
remote sensing and physiological ear and hearing systems to necessary parallel pro- 
cessing of large frame systems. Fusion frames provide a tool for weighted information 
combination from a set of overlapping subspaces in a distributed manner. The fusion 
process can be either based on local subspace processing/reconstructions through a 
fusion frame operator, or based on a global reconstruction with a distributed dual 
frame evaluation. More notably, for fusion frame systems, the fusion frame operator 
itself is numerically simple and efficient. It is a natural combination of local frame 
operators, which can be expressed as sums of local frame matrix representations. This 
makes fusion frames more than a mere notion but a practically and numerically handy 
tool for distributed processing. Stability is also a substantial feature of fusion frames, 
which makes it a robust distributed fusion processing system. 

For distributed systems such as general sensor networks, fusion frames are a tool 
easy and ready for system modeling and information combination. For large frame 
systems requiring parallel processing, fusion frames provide a means to subdivide the 
large system into a set of rather flexible and overlapping small subsystems. Each 
subsystem can be processed independently and then combined coherently. The com- 
putational efficiency is comprehensible. 

We envision that applications of fusion frames can reach afar with ample impact. 

Acknowledgments. Parts of the research for this paper was performed while 
the second author was visiting the Department of Mathematics at the University of 
Missouri in Columbia. This author thanks this department for its hospitality and 
support during this visit. First of all the authors are very thankful to J. C. Tremain 
for very helpful suggestions concerning the notion of perturbation for fusion frames 
which improved the paper. We are also indebted to H. Bolcskei, H. Feichtinger, J. 
Tanner, and K. Vasudevan for interesting discussions concerning further applications 
of fusion frames. 

REFERENCES 

[1] A. Aldroubi, C. Cabrelli, and U. M. Molter, Wavelets on irregular grids with arbitrary 
dilation matrices and frame atoms for L2(M<*), Appl. Comput. Harmon. Anal., 17 (2004), 
pp. 119-140. 

[2] A. Aldroubi and K. Grochenig, Beurling-Landau-type theorems for non-uniform sampling 
in shift invariant spline spaces, J. Fourier Anal. Appl., 6 (2000), pp. 93-103. 

[3] M. S. AsGARi AND A. Khosravi, Frames and bases of subspaces in Hilbert spaces, J. Math. 
Anal. Appl., 308 (2005), pp. 541-553. 

[4] J. Benedetto, A. Powell, and O. Yilmaz, Sigma-Delta quantization and finite frames, IEEE 
Trans. Inform. Th., to appear. 

[5] B. BODMANN, D. W. Krids, and V. I. Paulsen, Decoherence-insensitive quantum communi- 
cations by optimal C* -encoding, preprint (2006). 

[6] H. Bolcskei, F. Hlawatsch, and H. G. Feichtinger, Frame-theoretic analysis of oversampled 
filter banks, IEEE Trans. Signal Processing, 46 (1998), pp. 3256-3268. 

[7] E. J. Cande:s and D. L. Donoiio, New tight frames of curvelets and optimal representations of 
objects with piecewise Cf^ singularities. Comm. Pure and Appl. Math., 56 (2004), pp. 216— 
266. 

[8] P. G. Casazza, Custom building finite frames, in Wavelets, Frames and Operator Theory 
(College Park, MD, 2003), Contemp. Math. 345, Amer. Math. Soc, Providence, RI, 2004, 
pp. 61-86. 

[9] P. G. Casazza and O. Christensen, Perturbation of operators and applications to frame 
theory, J. Fourier Anal. Appl., 3 (1997), pp. 543-557. 



Fusion Frames and Distributed Processing 



21 



[10] P.G. Casazza and J. KOVACEVIC, Equal-norm tight frames with erasures, Adv. Comput. 

Math., 18 (2003), pp. 387-430. 
[11] P. G. Casazza and G. Kutyniok, Frames of subspaces, in Wavelets, Frames and Operator 

Theory (College Park, MD, 2003), Contemp. Math. 345, Amer. Math. Soc, Providence, 

RI, 2004, pp. 87-113. 

[12] P. G. Casazza and G. Kutyniok, A generalization of Gram-Schmidt orthogonalization gen- 
erating all Parseval frames. Adv. Comput. Math., to appear. 

[13] O. Christensen, An Introduction to Frames and Riesz Bases, Birkliauser, Boston, 2003. 

[14] M. S. Craig and R. L. Center, Geophone array formation and semblance evaluation. Geo- 
physics, 71 (2006), pp. 1-8. 

[15] I. Daubechies, Ten Lectures on Wavelets, SIAM, Philadelphia, 1992. 

[16] M. Fornasier, Quasi-orthogonal decompositions of structured frames, J. Math. Anal. AppL, 
289 (2004), pp. 180-199. 

[17] V. K. GOYAL, M. Vetterli, and N. T. Tilvo, Quantized overcornplete expansions in R'^ : 
Analysis, synthesis, and algorithms, IEEE Trans. Inform. Theory, 44 (1998), pp. 16-31. 

[18] K. H. Grochenig, Accelaration of the frame algorithm, IEEE Trans. Signal Processing, 41 
(1993), pp. 3331-3340. 

[19] R. W. Heath and A. J. Paulraj, Linear dispersion codes for MIMO systems based on frame 

theory, IEEE Trans. Signal Processing, 50 (2002), pp. 2429-2441. 
[20] S. S. Iyengar and R. R. Brooks, eds.. Distributed sensor networks. Chapman & Hall/CRC, 

Baton Rouge, 2005. 

[21] S. Li, On general frame decompositions, Numer. Funct. Anal. Optim., 16 (1995), pp. 1181—1191. 

[22] S. Mallat, a wavelet tour of signal processing, Academic Press, Inc., San Diego, CA, 1998. 

[23] C. J. RozELL AND D. H. Johnson, Analysis of noise reduction in redundant expansions under 
distributed processing requirements, in International Conference on Acoustics, Speech, and 
Signal Processing, Philadelphia, PA, 2005. 

[24] C. J. RozELL ,'\ND D. H. .Johnson, Analyzing the robustness of redundant population codes in 
sonsory and feature extraction systems, Neurocomputing, to appear. 

[25] C. J. RozELL, I. N. Goodman, and D. H. Johnson, Feature-based information processing with 
selective attention, in International Conference on Acoustics, Speech, and Signal Process- 
ing, Toulouse, France, 2006. 

[26] T. Strohmer and R. W. Heath, Jr., Grassmannian frames with applications to coding and 
communication, Appl. Comput. Harmon. Anal., 14 (2003), pp. 257—275. 

[27] W. Sun, G-frames and G-Riesz Bases, J. Math. Anal. Appl., to appear. 

[28] W. Sun, Stability of G-frames, preprint (2006). 

[29] A. TOSELLI and O. Widlund, Domain decomposition methods — algorithms and theory. 

Springer Series in Computational Mathematics 34, Springer- Verlag, Berlin, 2005. 
[30] E. Waltz and J. Llinas, Multisensor Data Fusion, Artech House, 1990. 



